Automated Boundary Conditions for Wind Tunnel Simulations

Computational fluid dynamic (CFD) simulations of models tested in wind tunnels require a high level of fidelity and accuracy particularly for the purposes of CFD validation efforts. Considerable effort is required to ensure the proper characterization of both the physical geometry of the wind tunnel and recreating the correct flow conditions inside the wind tunnel. The typical trial-and-error effort used for determining the boundary condition values for a particular tunnel configuration are time and computer resource intensive. This paper describes a method for calculating and updating the back pressure boundary condition in wind tunnel simulations by using a proportional-integral-derivative controller. The controller methodology and equations are discussed, and simulations using the controller to set a tunnel Mach number in the NASA Langley 14- by 22-Foot Subsonic Tunnel are demonstrated

Inflow/Outflow Boundary Conditions with Application to FUN3D

Several boundary conditions that allow subsonic and supersonic flow into and out of the computational domain are discussed. These boundary conditions are demonstrated in the FUN3D computational fluid dynamics (CFD) code which solves the three-dimensional Navier-Stokes equations on unstructured computational meshes. The boundary conditions are enforced through determination of the flux contribution at the boundary to the solution residual. The boundary conditions are implemented in an implicit form where the Jacobian contribution of the boundary condition is included and is exact. All of the flows are governed by the calorically perfect gas thermodynamic equations. Three problems are used to assess these boundary conditions. Solution residual convergence to machine zero precision occurred for all cases. The converged solution boundary state is compared with the requested boundary state for several levels of mesh densities. The boundary values converged to the requested boundary condition with approximately second-order accuracy for all of the cases

Assessment of an Explicit Algebraic Reynolds Stress Model

This study assesses an explicit algebraic Reynolds stress turbulence model in the in the three-dimensional Reynolds averaged Navier-Stokes (RANS) solver, ISAAC (Integrated Solution Algorithm for Arbitrary Con gurations). Additionally, it compares solutions for two select configurations between ISAAC and the RANS solver PAB3D. This study compares with either direct numerical simulation data, experimental data, or empirical models for several different geometries with compressible, separated, and high Reynolds number flows. In general, the turbulence model matched data or followed experimental trends well, and for the selected configurations, the computational results of ISAAC closely matched those of PAB3D using the same turbulence model

FUN3D and USM3D Analysis of the Propulsion Aerodynamic Workshop 2018 S-Duct Test Case

This work presents the results of Fun3D and USM3D analyses that were performed for the 4th AIAA Propulsion Aerodynamics Workshop (PAW). The PAW workshop is separated into three sections that focus on internal duct flows, nozzle flows and a special topic. This paper focuses on the internal duct flow section of PAW04 while an accompanying paper discusses the analyses performed for the nozzle portion. For the internal duct flow section, the PAW04 participants were provided with the two configurations consisting of an S-duct with and without aerodynamic interface plane (AIP) rake legs modeled. The participants were asked to perform a grid refinement study as well as a turbulence model study for the configuration with the rake legs. The analyses discussed here were performed on custom grids developed under the guidelines of the workshop. Additionally, the paper discusses the development and use of flow controllers for matching the desired flow characteristics. The results show that both solvers do well for predicting internal flow characteristics of the S-duct based on direct comparison with the experimental data. However, the CFD-to-CFD comparison proved to be more challenging due to the localized occurrence of supersonic flow near the rake legs when using the mass flow controller. A turbulence model study was performed to compare the two-equation SST model to the SA-QCR model. The results show that although the turbulence model does affect the solution, it makes a minimal impact on pressure recovery and inlet distortion intensity for this case. Suggestions for future workshops include gridding guidelines similar to those employed for the Drag Prediction Workshop series for the grid refinement study and a time accuracy study

Turbulent Output-Based Anisotropic Adaptation

Controlling discretization error is a remaining challenge for computational fluid dynamics simulation. Grid adaptation is applied to reduce estimated discretization error in drag or pressure integral output functions. To enable application to high O(10(exp 7)) Reynolds number turbulent flows, a hybrid approach is utilized that freezes the near-wall boundary layer grids and adapts the grid away from the no slip boundaries. The hybrid approach is not applicable to problems with under resolved initial boundary layer grids, but is a powerful technique for problems with important off-body anisotropic features. Supersonic nozzle plume, turbulent flat plate, and shock-boundary layer interaction examples are presented with comparisons to experimental measurements of pressure and velocity. Adapted grids are produced that resolve off-body features in locations that are not known a priori

Validation of a Node-Centered Wall Function Model for the Unstructured Flow Code FUN3D

In this paper, the implementation of two wall function models in the Reynolds averaged Navier-Stokes (RANS) computational uid dynamics (CFD) code FUN3D is described. FUN3D is a node centered method for solving the three-dimensional Navier-Stokes equations on unstructured computational grids. The first wall function model, based on the work of Knopp et al., is used in conjunction with the one-equation turbulence model of Spalart-Allmaras. The second wall function model, also based on the work of Knopp, is used in conjunction with the two-equation k-! turbulence model of Menter. The wall function models compute the wall momentum and energy flux, which are used to weakly enforce the wall velocity and pressure flux boundary conditions in the mean flow momentum and energy equations. These wall conditions are implemented in an implicit form where the contribution of the wall function model to the Jacobian are also included. The boundary conditions of the turbulence transport equations are enforced explicitly (strongly) on all solid boundaries. The use of the wall function models is demonstrated on four test cases: a at plate boundary layer, a subsonic di user, a 2D airfoil, and a 3D semi-span wing. Where possible, different near-wall viscous spacing tactics are examined. Iterative residual convergence was obtained in most cases. Solution results are compared with theoretical and experimental data for several variations of grid spacing. In general, very good comparisons with data were achieved

FUN3D and USM3D Analysis of the 4th AIAA Propulsion Aerodynamics Workshop Nozzle Test Case

This work presents the results of FUN3D and USM3D analyses that were performed for the 4th AIAA Propulsion Aerodynamics Workshop. The workshop was separated into three sections that focus on internal duct flows, nozzle flows, and a special topic of interest. This paper focuses on the nozzle section while an accompanying paper discusses the analyses performed for the internal duct flow section. For the nozzle flow section, the participants were provided with two configurations consisting of a rectangular convergent nozzle with and without an aft-deck. User-generated grids were developed under the guidelines of the workshop for the present analyses. The results show that both solvers compare favorably to the experimental results for the baseline nozzle with the largest differences observed at the lower NPR values. Additionally, both solvers showed favorable agreement with the experimental data for the pressures on the surface of the aft-deck. However, neither solver was able to match the jet flow further downstream for the nozzle with aft-deck configuration. A turbulence model study was conducted to compare the two-equation SST model, SA-QCR model, and two-equation k-kL model (FUN3D only). The results show that the SA-QCR turbulence model was unable to match the experimental results downstream of the nozzle. The k-kL model was shown to better match the experimental data compared to the SST model for most cases simulated using the FUN3D flow solver. Suggestions for future workshops include gridding guidelines similar ion Workshop series for the grid refinement study and a reduction in scope to allow for more detailed exploration of the individual problems. rticipation by reducing time requirement

Verification and Validation of the k-kL Turbulence Model in FUN3D and CFL3D Codes

The implementation of the k-kL turbulence model using multiple computational uid dy- namics (CFD) codes is reported herein. The k-kL model is a two-equation turbulence model based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's equation is well suited for term-by-term mod- eling and displays useful features compared to other two-equation models. An important di erence is that this formulation leads to the inclusion of higher-order velocity derivatives in the source terms of the scale equations. This can enhance the ability of the Reynolds- averaged Navier-Stokes (RANS) solvers to simulate unsteady ows. The present report documents the formulation of the model as implemented in the CFD codes Fun3D and CFL3D. Methodology, veri cation and validation examples are shown. Attached and sepa- rated ow cases are documented and compared with experimental data. The results show generally very good comparisons with canonical and experimental data, as well as matching results code-to-code. The results from this formulation are similar or better than results using the SST turbulence model

Nonlinear Dynamic Modeling of a Supersonic Commercial Transport Turbo-Machinery Propulsion System for Aero-Propulso-Servo-Elasticity Research

This paper covers the development of an integrated nonlinear dynamic model for a variable cycle turbofan engine, supersonic inlet, and convergent-divergent nozzle that can be integrated with an aeroelastic vehicle model to create an overall Aero-Propulso-Servo-Elastic (APSE) modeling tool. The primary focus of this study is to provide a means to capture relevant thrust dynamics of a full supersonic propulsion system by using relatively simple quasi-one dimensional computational fluid dynamics (CFD) methods that will allow for accurate control algorithm development and capture the key aspects of the thrust to feed into an APSE model. Previously, propulsion system component models have been developed and are used for this study of the fully integrated propulsion system. An overview of the methodology is presented for the modeling of each propulsion component, with a focus on its associated coupling for the overall model. To conduct APSE studies the de- scribed dynamic propulsion system model is integrated into a high fidelity CFD model of the full vehicle capable of conducting aero-elastic studies. Dynamic thrust analysis for the quasi-one dimensional dynamic propulsion system model is presented along with an initial three dimensional flow field model of the engine integrated into a supersonic commercial transport

Simulations of the NASA Langley 14- by 22-Foot Subsonic Tunnel for the Juncture Flow Experiment

NASAs Transformational Tools and Technologies Programs Juncture Flow experiment aims to provide data to improve Computational Fluid Dynamics (CFD) modeling in the juncture flow region. The experiment is planned to provide validation-quality data for CFD that focuses on the separation bubble near the wing-body juncture trailing edge region. Because wind tunnel tests associated with the Juncture Flow project have been designed for the purpose of CFD validation, considerable effort is going into modeling and simulating the wind tunnel. This is not only important because wind tunnel wall effects can play a role in integrated testing uncertainties, but also because the better the boundary conditions are known, the better CFD can accurately represent the experiment. This paper builds on the recent CFD efforts to model the NASA Langley 14- by 22-Foot Subsonic Tunnel. Current best practices in simulating wind tunnels are evaluated. The features of each method, as well as some of their pros and cons, are highlighted. Boundary conditions and modeling techniques currently used by CFD for empty-tunnel simulations are also described. Preliminary CFD studies associated with modeling the Juncture Flow model are summarized, with the intention to determine sensitivities of the flow near the wing-body juncture region of the model to a variety of modeling decisions